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PACS. 67.40.Vs - Vortices and turbulence. 

PACS. 67.57.-z - Superfiuid phase of liquid helium-3. 



Abstract. - In a recent experiment Finne et al. discovered an intrinsic condition for the onset 
of quantum turbulence in 3 He-B, that q — a/(l — a') < 1, where a and a' are mutual friction 
parameters. The authors argued that this condition corresponds to Kelvin waves which are 
marginally damped, so for q > 1 Kelvin waves cannot grow in amplitude and trigger vortex 
reconnections and turbulence. By analysing both axisymmetric and non-axisymmetric modes 
of oscillations of a rotating superfiuid, we confirm that in the long axial wavelength limit the 
simple condition q — 1 is indeed the crossover between damped and propagating Kelvin waves. 

o 

Introduction. - A striking feature of superfluidity is the existence of quantised vortex 
filaments. Quantised vortices are particularly interesting because, unlike vortices in a classi- 
cal turbulent fluid, they are discrete, stable topological defects. Until recently the study of 
quantised vorticity and quantum turbulence [1] was mainly confined to 4 He, but in the last 
, few years there has been progress in other contexts [2], ranging from atomic Bose-Einstein 
condensates to 3 He, which is a fermionic superfiuid. 

Our work is motivated by the recent discovery by Finne et al. [3] of a sharp transition at 
the temperature T — 0.6 T c between two distinct hydrodynamical regimes in 3 He-B, where 
T c is the critical temperature. In their experiment Finne et al. injected one or a few seeding 
vortex loops into rotating 3 He-B. Using a non-invasive NMR measurement technique, they 
measured the total length of the quantised vortices as a function of time. They found that 
if T > 0.6 T c the initial loops expand and straighten into rectilinear vortices which align 
along the axis of rotation. However, if T < 0.6T C , the seeding loops evolve into a turbulent 
vortex tangle. The effect is independent of the fluid's velocity and the details of the initial 
condition, which suggests the existence of a remarkable intrinsic criterion for the onset of 
quantum turbulence. On the contrary, the condition for the onset of turbulence in a classical 
fluid, that the Reynolds number is large, is extrinsic, because the Reynolds number depends 
on the geometry and the scale of the flow, not only on properties of the fluid. 

Finne et al. interpreted their experimental result in terms of Kelvin waves on the vortex 
filaments - helical perturbations of the position of a vortex core away from the unperturbed 
straight shape (see Fig. 1). They argued that the superfiuid vorticity is determined by the 
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ratio q of dissipative and inertial forces in the superfLuid, in analogy with the classical Reynolds 
number, which expresses the ratio of inertial and viscous forces in an ordinary Navier-Stokes 
fluid. The quantity q is defined as 

where a and a' are known [4] mutual friction coefficients, which determine the strength of the 
interaction between the superfluid vortices and the normal fluid. Finne et al. noticed that, for 
a single superfluid vortex, it was predicted [2] that q = 1 is the crossover from Kelvin waves 
which propagate (q < 1) and Kelvin waves which are overdamped (q > 1). By performing 
enlightening numerical simulations, Finne et al. confirmed that, if q > 1 an injected vortex 
loop straightens and becomes rectilinear, whereas if q < 1 Kelvin waves propagate freely, as 
predicted. In the latter case, the amplitude of the Kelvin wave grows, driven by the local 
difference between the normal fluid velocity and the velocity of the vortex line [5] . When the 
waves' amplitude becomes of the same order of the average intervortex spacing, the vortices 
reconnect with each other, quickly forming a disordered turbulent tangle. The value of q 
depends on T, and the observed temperature 0.6 T c corresponds to q = 1.3, which is indeed 
close to unity. 

The aim of our work is to reconsider the condition q = 1 (which for a single vortex is the 
crossover from free to overdamped motion) for a large density of vortices, such as the vortex 
lattice in the rotating helium experiment of Finne et al. . In this the context the behaviour 
of the superfluid is described by the Hall-Vinen equations [6]. Thus, rather than individual 
vortex lines, we consider a continuum of vortex lines. The advantage of this model is that 
it allows us to explore effects which the theory of a single vortex filament cannot describe, 
notably the presence of boundaries and the degrees of freedom represented by the coherent 
oscillatory motion of many vortices. Mathematically, this corresponds to the fact that for 
a single vortex, since the core radius is fixed, only m = 1 modes (sideway displacements of 
the vortex core) are possible, whereas for a bundle of vortices any azimuthal m symmetry is 
permitted. The simple question which we address is therefore the following: given a superfluid 
in a cylindrical container of radius a rotating at constant angular velocity f2, thus forming 
a vortex lattice with a large density of vortex lines, does the condition q = 1 represent the 
crossover from propagating Kelvin waves to overdamped Kelvin waves for all possible modes 
of oscillations ? 

The governing equations and the basic state. - The equation of motion of the superfluid 
in a coordinate system rotating with angular velocity $1 = fle z may be written as 

dv s 

— + (v s ■ V)v s = V* + 2v s x n + aX x [A x (v s - v n )] 

dt (2) 

+ a'A x (v s - v n ) - av s \ x (A • V)A + v s {l - a'){\ ■ V)A 

where v s and v n are the superfluid and normal fluid velocities in the rotating frame, A = 
Vxo s + 2S~i, A = A/|A| is the unit vector in the direction of A and \1/ is a collection of scalar 
terms. Given the high viscosity of 3 He-B, we assume that the normal fluid is in solid body 
rotation around the z-axis, thus in the rotating frame v n = 0. Eq. @ must be solved under 
the condition that V- V s = 0. The quantity v s = (T/Air) log(6 /a ) is the vortex tension 
parameter, T is the quantum of circulation, a is the vortex core radius and b = (|A|/r) -1 / 2 
is the average distance between vortices. The unperturbed vortex lattice corresponds to the 
basic state v$ = vfi = 0, V^o = for which A = 2fle z , working in cylindrical coordinates 
(r, <f), z). 
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The equations of the perturbations. - We perturb the basic state by letting v s = u = 
(u r ,u,p,u z ), \f r = + 4>, where \u\ <C 1 and \if}\ -C 1 are small perturbations. We assume 
normal modes of the form exp(icrt + imcf) + ikz)], where to and k are respectively the azimuthal 
and axial wavenumbers and a is the growth rate. The aim of our calculation is to determine the 
real and imaginary parts of er, namely TZe(a) and Tm(o~). The resulting linearised equations 
for the perturbations are 



/\ dip mkv . du z 

(la + ari)u r — 77(1 — a )iu = — (1 — a )u z — ikv a— p- = 0, (3) 

dr r dr 

, mi amkv s ,.du z 

T](l — a )u r + (i<7 + ocqjUA, — — ip H u z — \kv (1 — a. )—— = 0, (4) 

r r dr 

\gu z = lkip, (5) 

1 d du r 1m 

- — (r-j— )H u <j> + \ku z = 0, (6) 

r dr dr r 

where 77 = 2f2 + v s k 2 . Eliminating u r , and ip we obtain the following differential equation 
for u z : 

1 d ( du 

where 



r 

r dr V dr 



I3 2 u z = 0, (7) 



,2 = -fc 2 [(ia + ^) 2 + (l-aQV] 

' [{10- + ar])(io- + k 2 v s a) + {1 - a') 2 v s k 2 n} , [ 1 

The solution of Eq. (7J which is regular as r — > is the Bessel function of the first kind of 
order to, J m (j3r). To determine (3 we enforce the boundary condition u r — on the wall of 
the container r = a, which yields the secular equation 

{kaf JjJM [{la + av) 2 + a , )W] + 2mn(j{l al) = „ (Q) 
{(3a) J m (pa) 

where the prime denotes the derivative of the Bessel function with respect to its argument. 

Limiting cases. - The dispersion law which we have obtained, represented by Eq.'s I|8l9fl . 
has two limiting cases which have already appeared in the literature. Firstly, if we set a = 
a' = v s = 0, we recover Chandrasekhar's classical result for the modes of vibrations of a 
rotating column of liquid [7] : a is determined by 

/3 2 = k 2 (^L-l), (10) 

where the values of (3 are given by the roots of 

af3aJ' m {f3a) + 2mnj m (f3a) = 0. (11) 

Secondly, if we set a = a' = m = 0, Eq. (jSJ) reduces to J' {(3a) = 0. Let be the j th zero of 
JliO = J o(0 ( that is 6 = 3.83171, 6 = 7.01559 etc); then we have 
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(k 2 a 2 i] 2 + v s rjk 2 ^ 2 ) 
{& + k 2 a 2 ) ' 



(12) 
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If the wavelength of the perturbations is larger than the radius of the cylinder we may neglect 
edge effects, (ka) 2 ^> £?, and we recover Hall's result [8] for the axisymmetric oscillations of 
a rotating superfluid: 

a 2 w rj 2 = 2Q + v s k 2 . (13) 

Eq. l(T5j l generalises to a vortex lattice the dispersion relation a — v s k 2 of an individual vortex 
line [2]. 

The axisymmetric case. - Eq.'s 18191 simplify in the axisymmetric case, m = 0. We have 
again Jg(/3a) = for which /3a = £j and Eq. becomes 



t] [(i<7 + my) 2 + (1 - a') V 



(fca) 2 [(id + a77)(i(j + av s k 2 ) + (1 — a') 2 v s rjk 2 ] 
This quadratic equation has solutions 

((afc) 2 + g)a = ia(2n(ak) 2 + tig + v s k 2 ({akf + rf)) 



(14) 



.(15) 



± ^/{(l - a') 2 k 2 {(ak) 2 + S, 2 )(2fl + v s k 2 )(2fla 2 + v s {{ak) 2 + £ 2 ))) - a 2 ft 2 ^ 4 } 

This result is identical to that derived by Glaberson et al. [9] who expanded in normal modes 
of the form exp(ifcix + ik%y + lkz + lat) where /a 2 = kf + k\. Expanding in normal modes 
of this form means that the boundary condition at r — a is not enforced. 

From Eq. I|15[) it can be seen that Tm(a) is non-negative, so the system is always stable to 
infinitesimal disturbances. If we neglect perturbations in the non-axial direction (setting = 
£ = fci = k% = 0) we obtain 

a = ia(20 + v s k 2 ) ± (1 - a')(2ft + v s k 2 ), (16) 

as quoted by Hall [8]. This solution is illustrated in Fig. [21 by a dashed line and satisfies 
\Tm(a) /lZe(a)\ = q. However, if we do not neglect perturbations in the non-axial direction, 
Eq. (|15fl has an infinite number of solutions, according to the value of considered. The least 
stable mode is the one for which Xm(a) is minimum and is the mode which we would expect to 
be observed providing the initial state is sufficiently random. Tm(a) decreases monotonically 
as k — > and £j — > oo. For small k the term inside the square root of Eq. (|15|l is negative, 
resulting in IZe(a) being zero and the mode decays monotonically with time. As k increases, 
the term inside the square root is positive and the mode displays oscillatory decay with time. 
The solution which has the smallest decay rate is 



a = ia(n + v s k 2 ) ± - a') 2 v s k 2 {2tt + v s k 2 ) - a 2 ^ 2 } (17) 

which corresponds to £j — > oo. Plots of the decay rate, —Xm{cr) and the ratio of the decay 
rate of the wave and its angular frequency, \2m(a) /lZe(a)\ can be found in Fig. [21 where the 
mutual friction parameters correspond to T = 0.4T C , for which a = 0.1125 and a' = 0.1, and 
we have taken v s /a 2 — 0.001 and O = 1. For large (ak) we find 

\lm(a)/ne(a)\^-^—=q. (18) 
(1 - a') 

The minimum permitted value of k will be governed by the height, h of the apparatus by the 
relation fc m i n = 2n/h, so provided the aspect ratio h/a is small enough we find agreement 
with the argument of Finne et al. [3] . 
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Non-axisymmetric case. In order to consider the non-axisymmetric modes we must 
solve the coupled equations I|8I9|I . For given temperature, T, wavenumbers m, ak and param- 
eters v s jo? , f2 Eq.'s (18191 were solved using NAG routine C05NDF on the real and imaginary 
parts of a and (a/3) 2 . As for the axisymmetric case, there will be an infinite number of 
solutions for each mode, m considered. 

In Fig. we illustrate some of the solutions computed for the m = 1 mode. The quantities 
— Zm(er) and \Tm(a) /lZe(a)\ are plotted against ak using the same parameters as for Fig.|3 
The dotted line represents the least stable axisymmetric mode. We can see that the least 
stable m = 1 mode is bound by the least stable axisymmetric mode and that, as for the 
m = case, all the computed solutions show \Tm(o~) /lZe(o~)\ — > q at large values of ak. 

In Fig. Q] we plot — lm(a) and \Tm(a)/TZe(a)\ against ak of the least stable computed 
modes for m = 0,2,4,8, 11, oo using the same parameters as for Fig. [21 The arrow shows 
the direction of increasing m. From the values of Im(<j) it can be seen that all the modes 
are stable and that the least stable mode occurs as m — > oo. For all the modes we find that 
\Im{a) /lZe(a)\ — > q for large values of ak and using the same argument put forward for the 
axisymmetric case, we conclude that we find agreement with the argument of Finne provided 
the aspect ratio of the apparatus is small enough. 

As m — > oo we are able to make a simplification to Eq. @. From Abromavitch & Ste- 
gun [10] we find 

1 / ez \ m 

J m {z) « i= hq— as m ^ oo, (19) 

V27tto \2m/ 



so, in this limit, we may take 



(a(3)J m (aP) " JaW' 

Substituting into Eq. © we find that the resulting equation is independent of m. This 
equation may be solved analytically and yields two solutions for a, namely a = (1 — a')v s k 2 + 
iav s k 2 and a = —(1 — a')r) + iarj. The first solution will be the least stable and this is 
the one that has been plotted in Fig. Note that both of these solutions are such that 
\lm(a)/Ke(a)\ = q. 

Conclusion. - In conclusion, our analysis shows that in the case of both axisymmetric 
(m = 0) and non-axisymmetric (m ^ 0) perturbations of a vortex lattice, provided that the 
aspect ratio is small enough (h/a ^> 1), Kelvin waves propagate if q < 1 and are damped if 
q > 1. This result confirms the argument put forward by Finne et al. . Finally, the exact 
dispersion relation which we have found (Eq. |Sland|^J) can be used to study with precision the 
spectrum of Kelvin waves on a rotating vortex lattice for any height and radius of experimental 
apparatus. 
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Fig. 1 - (a): unperturbed vortex; (b): Kelvin wave. 



(a) (b) 




Fig. 2 - Plots of (a) the decay rate, Tm(a) (b) \lm(a) /TZe(a)\ of the axisymmctric mode (m = 0) 
against ak for various values of The arrows show the direction of increasing where j = 
1, 2, 5, 10, 20, 100, oo. The dashed line shows the result for £j — 0. 



(a) (b) 




Fig. 3 - Plots of (a) the decay rate, Im(u) (b) \Tm(a) /lZe(a)\ against ak of the m = 1 mode. The 
dotted line represents the least stable axisymmetric mode. 
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Fig. 4 - Plots of (a) the decay rate, Tm(a) (b) \Im(a) /lZe(a) \ against ak for the m = 2, 4, 8, 11 mode. 
The dotted line represents the least stable axisymmetric mode. The solid line represents the m — oo 
mode. The arrow shows the direction of increasing m. 
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